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ABSTRACT 

The ultrahigh-energy cosmic-ray anisotropics discovered by the Pierre Auger Observatory give the potential 
to finally address both the particles' origins and properties of the nearby extragalactic magnetic field (EGMF). 
We examine the implications of the excess of ~ 10^^ eV events around the nearby radio galaxy Centaurus A. 
We find that, if Cen A is the source of these cosmic rays, the angular distribution of events constrains the 
EGMF strength within several Mpc of the Milky Way to > 20 nG for an assumed primary proton composition. 
Our conclusions suggest that either the observed excess is a statistical anomaly or the local EGMF is stronger 
than conventionally thought. We discuss several implications, including UHECR scattering from more distant 
sources, time delays from transient sources, and the possibility of using magnetic lensing signatures to attain 
tighter constraints. 
Subject headings: ISM: cosmic rays — ISM: magnetic fields — galaxies: magnetic fields 



1. INTRODUCTION 



Significant progress has been made toward determin ing the 
origin of ultrahigh-energy cosmic rays (UHECRs; see 
T9Mand 'Beatty & Westerhoff20Q9^, 'Kotera & Olinto 



Hillas 



2011 



Letess ier-Selv on & Stan ev 201 1 for recent reviews), with the 
Pierre Auger Observatory (Auger) acquiring data needed to 
achieve this and other goals. However, a definitive identifica- 
tion of sources remains elusive, as early indications of a cor- 
relation of UHECR events with AGN within ^100 Mpc from 



the Veron-Cetty and Veron catalog (Pierre Auger CoUabora- 
|tion 2007) have given way t o a less-clear, y et still anisotropic, 
picture (Pierre Auger Colla boration|20 1 Oal) . 

The most prominent feature in the UHECR sky is a sig- 
nificant excess of even ts from the vicinity of Centaurus A 
(Cen A; see Fi g. [T]) (Pier re Auger Collaboratio n 2008, 2010a; 
Gorbunov et al. 2008; Stanev 2008; Moskalen ko et al.|2009t 
Hillas,2009) , a nearby acti ve radio galaxy pos s essing well- 
studied gian t radio lobes ( [Junkes et al.||1993] |Feain et aL 



2009 1 12011"). These radio lobes, along with the proximity to 



the Milky Way, ha ve long made Cen A a prime prospective 
source of UHECR ([Cav allo 1978 ; Romero et al.n996t [Ahn| 
^FaL]|T999l IFarrar & Piran. 2000, .Isola et al.. 2002i rThis is 
in stark contrast to the lack of an excess towards the Virgo 
galaxy cluster, home of the powerful AGN M87. 

A major difficulty in tracing cosmic rays, even those with 
energies exceeding 10^^ eV, back to their sources is the uncer- 
tain n ature of both the Galactic magnetic field (GMF) (Stanev] 
1991) and the extragalactic magnetic field (EGMF) in th e 



vicinity of the Milky Way ( lKronberg|p^94| |Widrow|[2Q02| ). 
The EGMF remains poorly constrained, especially locally. 
While upper limits of ~ 1-10 nG with coherence lengths of 

^ Theoretical Division, MS B285, Los Alamos National Laboratory, Los 
Alamos, NM 87545 

^ Bartol Research Institute, Department of Physics and Astronomy, Uni- 
versity of Delaware, Newark, DE 19716 

^ Lawrence Berkeley National Laboratory and Department of Physics, 
University of California, Berkeley, CA 94720 

4 California Institute of Technology, MC 350-17, Pasadena, CA 91 125 

^ Einstein Fellow 

^ Department of Physics, University of Toronto, Toronto M5S 1A7, 
Canada 



^ 1-50 Mpc have been placed in cosmological settings CB lasi| 
|et al.|19 99), a global estimate of the energy output of all su- 
permassive black holes implies '^0.1-1 jiG fields sp read out 
within the fi laments connecting clusters of galaxies ( jColgate| 
et al.||2011 ). Faraday rotation studies sugges t fields as large 
as - 0.3 fiG in the filaments (Xu et al.'2006') and - 1-10 /iG 
in intra-cluster environments (Clarke et al. 2001). A hope of 
UHECR studies is to increase our knowledge of the proper- 
ties of t he EGM F (e.g., .Kronbergi 19941 |Lee et al.|1995t[SigI 
etal.|19 99; Sta nev et al.|2000||Dolag et al |2005[|Anchordo- 
qui & G oldberg 2002VAnchordoqui et al."2001VArmengaud 
et al.||2 005; Ryu et al. 2010; Jiang et al. 2010; Anchordoqui 
et al. 2011 ; Takami et al. 2012). One way to do so requires 
evidence that a known object is indeed producing UHECRs in 
order to reconstruct the intervening field structure. 

We examine the implications of the Cen A excess by com- 
paring the observed distribution of > 5.5 x 10^^ eV Auger 
events with the expectations from UHECR propagation in a 
turbulent EGMF. Using this technique, we find that the over- 
all angular distribution can be well reproduced for a range of 




Fig. 1. — The arrival directions of 69 UHECR events detected by Auger 
with E > 5 .5 X IQ-*^^ eV (black circles) in Galactic coordinates ( Pierre] 
[Auger Collabo ration 2010aL Event pairs within 5° of each other are fur- 
ther crrcled. The light solid line demarcates the horizon of Auger, while the 
darker solid line corresponds to the supergalactic plane. We also show an es- 
timated density distribution (contours) obtained by replacing discrete events 
with a Gaussian and weighting by the relative exposure of Auger. A circle of 
18° radius is drawn around the center of the radio galaxy Cen A. 
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magnetic field configurations if Cen A is an UHECR source. 
This then allows us to constrain the strength and the structure 
of the local EGMF. Our focus is on the effects of the nearby 
EGMF on the propagation of UHECRs, and later we address 
why the GMF should not impact our conclusions. 

Throughout this study, we assume that Cen A is the source 
of the excess, and that the UHECR are protons. We will 



re-we ight each event using the exposure function ( Sommers 
200 1| ) at its location in the sky, so that events with lower ex 



discuss alternatives, while addressing Auger results ( Pierre 
Auger Collaboration] |20 1 1 1 ) that do not display the excess 



around Cen A that is expected at lower energi es, if the high- 
energy events are heavy nuclei ( Lemoine & Waxman|[T 009). 
The observed rapid isotropization with decreasing energy is 
most easily interpreted as proton domination and arises in 
subtle ways in a number of our viable field configurations. We 
develop further immediate implications of our scenario, in- 
cluding the effect on angular distributions from more distant 
sources, a requisite minimum time delay for UHECRs from 
transient sources (e.g., gamma-ray bursts), and the possibil- 
ity of using magnetic lensing effects to narrow the allowed 
magnetic field parameter space. 

2. CEN A'S PLACE IN THE UHECR SKY 

The deflection of UHECRs during propagation, which de- 
pends on their charge/momentum and any magnetic fields 
encountered, has long frustrated attempts at associating ob- 
served events with well-known astronomical objects. Ow- 
ing to their extreme energies, > 10^^ eV cosmic-ray protons 
should not experience strong deflections in typical Galactic 
(~ /iG) magnetic fields, and may point back to near their 
birthplaces if extragalactic fields are not too large ( Stanev et 
|al. 2000). A corollary to this is that the UHECR distribution 
from a known source can be used to determine the structure 
and the strength of the magnetic fields. 

We show in Fig. [T] the 69 UHECR events wi th £^ > 5.5 x_ 
10^^ eV (dots) reported through March 2009 (Pierre Auger 



Collaboration^ 2010a). In order to better visually represent 



the UHECR density distribu tion, we hav e also replaced each 
event with a Gaussian ( Silverman! 1986| ) of width 10° (com- 
parable to the spread of excess events around Cen A). To cor- 
rect for the non-uniform exposure of Auger, we also inversely 
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Fig. 2. — The cumulative angular distribution of the 69 Auger UHECR 
events when centered on Cen A (steps). We compare, after weighting for ex- 
posure, the expectations for: (i) a purely isotropic distribution of all events 
(solid), and; (ii) 10 events following a 10° Gaussian distribution around 
Cen A upon an otherwise isotropically distributed 59 events (dotted). 



posure display higher significance. The resulting density dis- 
tribution is normalized to unity and is displayed as shaded 
contours. 

It is clearly apparent that no other strong excess exists other 
than in the direction of Cen A. Quantitatively, 13 events are 
observed within 18° of Cen A, while ~3 events are expected 
on average for a pur ely isotropic distribution ( Pierre Auger| 
Collaboratio n" "2010a), suggesting an excess of about ^10 
events. In Fig. [2j we show the cumulative distribution of the 
Auger events in terms of the angular separation from Cen A 
(using Cos OcenA for uniform coverage in solid angle), along 
with the expectations, after weighting with the Auger expo- 
sure, for: (i) an isotropically distributed set of events (solid 
line); (ii) a 10° Gaussian containing 10 events centered at 
Cen A superimposed on an otherwise isotropic distribution 
(dotte d line). 

The [Pierre Auger Collaboration! ( |2010a| determined that 
only 4% of random realizations of 69 events drawn from an 
isotropic distribution deviate from the assumption of isotropy 
by more than the data itself do. Following this prescription, 
we find ^6% deviations for model (i), in good agreement with 
the [Pierre Auger Collaboration] ( |2010a |, and ~8% for model 
(ii). This suggests that even accounting for an excess from 
the direction of Cen A does not alter the fact that the all-sky 
distribution of events is anisotropic. Indeed, examining Fig.[T] 
we see a large low-density void, where the exposure is nearly 
maximal, that rivals the Cen A excess in prominence. 

The projected angular distribution is part of a larger story. 
Closer examination reveals that many of these events are con- 
tained in clusters, and further that most of the clusters are near 
Cen A. In Fig.[T] we put small circles around the events within 
5° of each other. Additionally, the events are not symmetric 
about the position of Cen A. These peculiarities may result 
from the particular structure of the EGMF, a possibility that 
we later address. 

There are a number of possibilities for accelerating UHECR 



in a radio galaxy like Cen A (Colgate 2004 ; Hardca stle et al. 
20091 IFraschetti & Mefia ;'2Q08 ; O' Sullivan et aL||2009[ 



mer et al.|| 2009| JRieger & A haronian 2009; Pe'er & 



Der- 



.oeb 



2012; Sahu et al.|2012| ). Acceleration in the inner jets or near 
the black hole both would occur on a scale much smaller than 
the radio lobes. It is also possible that particle acceleration 
occurs within the giant lobes themselves. For our purposes, 
we assume that for all instances particle injection arises from 
a central point source, since the apparent extent of the excess 
exceeds even the '^ 8° x 2° angular size of the lobes. 



3. MAGNETIC FIELD FORMULATION 

Any treatment of cosmic ray propagation in the nearby uni- 
verse requires a description of the magnetic field structure of 
the local EGMF. To simulate this, we generate a divergence- 
free, random magnetic field model whose components have 
a Gaussian distribution and follow \Bk\ ex k~^^~^'^\ Here 
n = 5/3 corresponds to Kolmogorov turbulenc e. Adopting 
the prescription detailed in |Tribble| ( |19 91) and Mur gia et al. 
(2004), we first construct Fourier components of a complex 
valued vector potential A whose components follow a power 
spectrum of the form \Ak\ ex k~^^~^^^ in a 3D cubical box in 
wave number (k) space. For a given k, each component Ai is 
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drawn from the distribution 

A 
P{Ai, (j)) dAi dcj) = — ^ exp 



-2 



B, [nG] 



A? 



2\A. 



dAi 



d4 

27r' 



(1) 



Solving for Ai, when the power spectrum normahzation is 
omitted, yields 



y^A;-("+4)ln[l//j2]e*2^p2^ 



(2) 



in which pi and p2 are real random numbers in the range 
to 1 . Construction of each A vector thus requires six random 
numbers. The values of A are calculated for all k in the range 

/Cmin = 27r/Cax < |^| < 27r /Cin = /^max, in which ^min 

and £max correspond to the minimum and maximum scales of 
the magnetic field. The magnetic field in k-space is given as 



B(k) = ikx A(k), 



(3) 



and is transformed back into coordinate space through a 3D 
complex fast Fourier transform. 

Fig.[3]shows one slice from the magnetic field configuration 
obtained using a cubic grid of size 512^. The z component of 
the field is projected in the x-y plane. Here, we choose /Cmin = 
3, /cmax = 256 for a Kolmogorov spectrum (although A for 
^max > 128 is assumed to be as the size of the grid will be 
too coarse to properly resolve them and these very small scale 
fluctuations would be under sampled). The quadratic mean 
of the magnetic field is normalized to ^rms = 1 nG over the 
volume and the wave numbers are scaled such that ^max = 
2 Mpc (for fc = 3) and ^min ^ 0.04 Mpc (for k = 128) in 
coordinate space. 

The most relevant quantities to describe any configuration 
of the magnetic field struct ure are Brms and the coherence 
length Ac, which is defined ( Harari et al.|2Q02| ) as 



A. 



L 



i i (,^min/^max/ 



77/ i V^min/^maxJ 



-1 • 



(4) 



When we use this magnetic field configuration to track 
the trajectories of charged particles for a given field strength 
(^rms) and coherence length (Ac), we renormalize B^ms and 
rescale ^min and ^max in our magnetic field simulation box ac- 
cordingly. Choosing kmin = 3 ensures that there is sufficient 
variation in the simulation box even at the largest scales. We 
take this field to be static over the relatively short propaga- 
tion times of the particles. The simulation box is periodically 
repeated when the propagation distance exceeds the box size. 

4. OBTAINING ANGULAR DISTRIBUTIONS 

In order to understand the angular distribution of events 
seen by Auger, we first look for a range of EGMF parame- 
ters that can result in a spread of ~ 10° for UHECR arriving 
from Cen A. We will use analytical expressions for guidance 
and later compare to in-depth simulations. In our equations, 
unless stated otherwise, we express both distance and time in 
Mpc (since v o:^ c for ultra-relativistic particles), energy in 
EeV (10^^ eV), magnetic field strength in nG, and consider 
singly-charged particles, applicable for protons. The equation 
of motion is governed by the Lorentz force, which, in a re- 
gion lacking electric fields, can be cast as a set of first order 
coupled equations 
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Fig. 3. — A slice from the magnetic field configuration, obtained from a 
Kolmogorov spectrum within a 512*^ cubic grid, used in this study. Shown is 
the z component of the field projected onto the x-y plane. 

where we express the velocity vector in terms of c, thus, (3 is 
a unit vector following these equations. The Larmor radius in 
a constant magnetic field is given as 



TL - 1.08 E/B. 



(6) 



While the energy of the particle in these equations can be re- 
placed with the rigidity R = E / Z for a. nucleus with charge 
Z, we keep this implicit and simply use E rather than R un- 
less otherwise noted. All of our subsequent results are quoted 
for a given energy and can be scaled from a proton to an iron 
primary by accounting for the appropriate charge through ei- 
ther the normalization of magnetic field or energy, although 
care must be taken in reinterpreting our ultimate conclusions 
for the case of heavy nuclei, as we discuss later. Due to the 
proximity of Cen A, we do not include energy losses here. 

In a turbulent magnetic field with a given strength, 5rms. 
and coherence length, Ac, the quadratic mean of the scattering 
of the final particle velocity with respect to the initial velocity 
is 

S,^, ^ 53° v/l72 ^rms VdA~c IE , (7) 



where d is the distance travelled ( [Harari et al.||2002| ). The 
corresponding scatter in arrival directions of particles around 
the source as seen at Earth is 

^rms = ^rms/v3. (8) 

While Eq. ^ is valid for d :^ A^ when d <C Ac the parti- 
cle effectively travels in a domain in which the magnetic field 
is constant. In that case, the deflection could be larger and 
will appear as a shift from the position of the source rather 
than scattered around the source. In this case, the change in 
the direction of the particle velocity, when averaged over all 
magnetic field configurations, can be expressed as 



(9) 



53° V2/3 5rms^/^, 
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Fig. 4. — The average angular distribution of 60 EeV cosmic-ray protons 
around Cen A. Shown are the expectations from Eq. flT) {dotted lines) and 
the results of simulations through turbulent magnetic fields {solid lines). 
We identify regions where the behavior is dominated by either scattering 
(Ac <C d) or shifting (Ac ^ d) of the source flux as seen by an observer. 
Magnetic lensing effects are strongest along the gray band (see text). 
Detailed results of simulations are presented for parameters denoted by large 
orange dots in Fig. [6] 



in which the factor y^2/3 accounts for the average perpendic- 
ular component of the random magnetic field direction with 
respect to the velocity. In this case, the shift in the observed 
arrival directions around the source is 



One can combine these two relations as 



77/2. 



1 



l/ry 



(10) 



(11) 



which describes either the average shift or the quadratic mean 
angular scattering of particles from the source. Here r] ^ —A 
parametrizes the smoothness of the transition between the two 
regions in which Kc^ d and Ac <C d. 

We consider particles with £^ = 60 EeV, comparable to the 
majority of the Au ger events, propag ating from a distance to 
Cen A of 3.8 Mpc ( Harris etal.|201Q| ). Fig.|4]shows values for 
(dotted lines) based on Eq. \\\\ for a given Ac and ^rms- We 
next compare this analytic result with numerical simulations. 

We compute 9 utilizing a fourth-order Runga-Kutta method 
to solve Eq. ([5]), keeping the step size small in comparison to 
both the minimum scale of variation in magnetic field and tl • 
In order to ensure that our results are always an average of 
many distinct realizations of the magnetic field configuration, 
we have varied the relative locations of the source and detec- 
tor while keeping their geometrical arrangements fixed when 
necessary. 

As seen in Fig.|4j the agreement between numerical results 
(solid lines) with the analytical results (dotted lines) is excel- 



lent. However, these results only provide the average of the 
distribution, losing important details from the azimuthal dis- 
tribution for comparison with the Auger data. Fig. [4] will be 
used as a guide as we pursue this direction via simulation. 

5. SIMULATIONS & SCENARIOS 

One can identify various realms, two of which have already 
been mentioned and labelled in Fig. |4] Scattering is charac- 
terized by Eqs. ([7]) & ([5]) and corresponds to a random walk 
around the source position for small angles. Larger scattering 
optical depth eventually leads to diffusive behavior. Shifting 
describes the case when the coherence length of the field ex- 
ceeds the distance to the source, so that the particles see an ap- 
proximately uniform field. This results in an apparent change 
in the position of the source as described by Eqs. ^ & ([TO]). 

Lensing refers to the magnification and/or appearance of 
multiple images of the cosmic-ray source due to uncorrelated 
deflections in the magnetic field. It roughly corresponds to the 
transition between the two realms mentioned above. This is 
an attractive possibility for using the observed angular distri- 
bution to infer greater details about the intervening magnetic 
field. For distant sources, the lensing effects are strongest near 
the critical energy. 



(12) 



Ec- 0.6B,^,d^^y^/A, 



as given in'Hara ri et al.| ( |20Q2] ). The gray band in Fig.|4]corre- 
sponds to the values of Brms versus Ac that produce the most 
prominent lensing from a source at a distance ofd^ 3.8 Mpc 
for particles of E c^ Ec/V^ ^ 60 EeV. Below this band, 
shifting dominates and we only have a single image, while 
above scattering dominates. Around this band, the source 
could be magnified, which may contribute to the excess of 
events from the direction of Cen A. 

To better demonstrate these realms, we have chosen mag- 
netic field parameters for UHECR propagation simulations in 
order to display realizations in which the all-sky averaged dis- 
placement of all trajectories from the center of Cen A is ~ 10° 
(denoted by the large orange dots in Fig. [4]). The top panel of 
Fig. [5] shows this for one such set. In these simulations, we 
have used a small test sphere at a distance of 3.8 Mpc from the 
source on which incoming particles are detected and we make 
the resulting sky plots. The average scatter (~ 2°) among the 
dark blue points in the lower right realization of Fig. [6] illus- 
trates the effective resolution determined by the size of the 
sphere, chosen as a compromise between higher statistics and 
image quality. 

We display the results of these simulations in Fig. |6] In 
the top panel of each of the four realizations, we show the 
all-sky distribution of UHECR (the positions of particles after 
reaching a sphere of radius 3.8 Mpc) as seen by an observer 
located at Cen A for cosmic rays of energies 60 EeV (dark 
blue points), to compare with the Auger data, and 10 EeV 
(light orange points) in order to illustrate behavior at lower 
energies. In the bottom panels, we show the sky distributions 
seen at Earth from three different locations relative to Cen A 
(as marked in the upper panels by arrows) for each simula- 
tion, along with the Auger data (black circles) for comparison. 
The location of Cen A is at the center of the large green cir- 
cle. These illustrate the unique features that can arise in these 
configurations. Despite this, the average angular spreads are 
roughly equal. 

Starting from small coherence lengths, we see in the up- 



Cen A UHECR Excess and Local Ext. Magnetic Field 



per panels that the distribution of trajectories as seen from the 
source is rather isotropic for the 60 EeV cosmic rays, due to 
the scattering being the largest. As Ac is increased, large-scale 
clustering becomes more prominent as the lensing regime is 
traversed. This gives way to a more isotropic dispersal again 
at large Ac due to the lack of scattering. For the 10 EeV cos- 
mic rays, a similar track is evident, although, since the EGMF 
parameters were chosen based on the 60 EeV results, this is 
not as well defined in the figures. For both energies, when 
lensing features are most prominent, significant fractions of 
the sky receive relatively low fluxes, so that a simple inverse- 
square based inversion may not yield the true luminosity of 
the source. 

As seen from Earth, in many cases when Ac is low a spread 
is present in the 60 EeV arrival directions, with sub-clustering 
among these sometimes present, along with an overall dis- 
placement from the center of Cen A. Conversely, larger co- 
herence lengths yield tight clustering with an overall shift in 
position relative to Cen A. These display how the average an- 
gle of all events in the simulation can be ~ 10°, as also given 
by the previous formulae, do not yet yield a distribution simi- 
lar to the Auger data around Cen A. This permits general con- 
clusions to be drawn regarding the properties of the EGMF 
despite the inherent difficulties involved with examining an 
infinite parameter space. 

6. THE LOCAL EXTRAGALACTIC MAGNETIC FIELD 

We have assumed in the above that the EGMF is the dom- 
inant source of the scattering of UHECRs from a nearby ob- 
ject and now provide constraints based on observations. We 
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Fig. 5. — Top: The observed angular distribution of 60 EeV cosmic-ray 
protons for EGMF parameters Brms = 16 nG and Ac = 1 Mpc, which 
yields an average displacement of ~ 10° from Cen A. Bottom: The corre- 
sponding distribution of delay times (propagation time minus d/c), averaging 
~ 0. 1 Myr. All EGMF parameters that are compatible with the UHECR dis- 
tribution around Cen A give similar results. 



address the deflections by GMF at the close of this Section. 
We show in Fig.|7]the regions in the parameter space of 5rms 
and Ac that satisfy two conditions based on the Auger data. 
First, we require the average angular distribution of events, 
averaged over the whole sky in order to use the highest statis- 
tics, to be 8-18° from Cen A. This yields the shaded region 
between the solid lines. 

As discussed above, large Ac values can meet this stan- 
dard yet still result in point-like distributions. We thus fur- 
ther require the "internal" spread among the arrival direc- 
tions to be larger than 4°, with our simulations resulting in 
the dashed lines. The intersection of these regions results in 
the preferred combination of EGMF parameters (dark shaded 
region). While for random field configurations, even with 
the same choice of parameters, it is difficult to define sharp 
boundaries on the allowed properties of the EGMF, we have 
attempted to estimate these under the assumption that Cen A 
is the cause of the UHECR excess coincident with its posi- 
tion in the sky. Techniques have been developed to compare 
in detail two-dimensional distributions of data (e.g.. Pea cock 
1983^), although, given the present data, our aim here is to 
provide a broad perspective. 

The most direct implication of this result is that, if Cen A 
is the source of protons resulting in the excess seen by Auger, 
then the strength of the intervening EGMF is > 20 nG. Mea- 
surements of the depth of maximum for UHECR showers at 
< 40 EeV made by Auger over large regions of the sky (i.e., 
not speci fically around Cen A), suggest a hea vy nuclei com- 
position ( [Pierre Auger Collaboration 2010b| ). Similar mea- 
surements from HiRes (located in the northern hemisphere 
and not covering Ce n A) instead indicate a light composition 
( Abbasi et al.||2010 ), which leads to an ambiguous situation 
concerning the all- sky composition. 

However, if the high-energy excess surrounding Cen A is 
due to heavy nuclei, this would imply a flux of protons of the 
same rigidity that propagate along the same trajectories if the 
accel erated particles were drawn from a solar-like composi- 
tion ( [Lemoine & Waxman||2009| ). This leads to an expected 
excess at lower energies of the same angular extent, though 
not as prominent due to the larger expected isotropic back- 
ground. This was not seen by Auger; however, leading to 
limits on the elemental composition to ratios that would be 
far fro m solar values if the > 55 EeV ex cess is due to heavy 
nuclei ( [Pierre Auger Collaboration|2011 ). The simplest inter- 
pretation of this result is a dominant proton component. 

This result has many consequences. A number of studies 
have suggested that the filling factor of extragalactic space 
containing fields of at least this strength is not very large (e.g., 
Sigl et al. 2004 ; [Dolag et al.[[2QQ5l [Das et al.[[2QQ8l [Takami 



et al.[ [2009i Gi acinti et al.[[201Q| ). However, the Milky Way 



being a fairly large galaxy may give credence to our being lo- 
cated within a filament containing a re latively- strong field, a s 
seen elsewhere in the nearby universe d Kronberg et al.|2007] ). 
Large scale simulations of extragalactic magnetic fields have 
difficulty in dealing with details at the scales of relevance 
here, which nonetheless can have large effects on UHECR 
propagation. 

From the available data alone, it is not possible to infer field 
properties beyond our neighborhood. However, a > 20 nG 
field extending at least > Mpc around the Milky Way results 
in a "screen" scattering all UHECRs that eventually reach 
Earth. Each UHECR would then be expected to have a mini- 
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Fig. 6. — Distributions of cosmic rays for four different extragalactic magnetic field parametrizations (Brms vs. Ac) using the magnetic field structure from 
Fig.js] Populations of protons with energies of 60 EeV (dark blue points), similar to the observed Auger event energies, and 10 EeV (light orange points) are 
shown. For each of the four models: Top: The final positions of particles, as seen by an observer at Cen A, after reaching a distance of 3.8 Mpc. 10^ particles are 
shown for each energy. Bottom: Three characteristic realizations of UHECR angular distributions arriving from Cen A, as seen at Earth, chosen from locations 
in the map above (as marked with arrows) shown together with Auger data (black circles). 



mum amount of deflection due to this field alone. If this were 
comparable to the angular extent of the Cen A excess, it may 
increase the difficultly of making associations with more dis- 
tant sources. Further, this local screen would introduce a min- 
imum time dispersion. To illustrate this effect, we show in 
Fig.|5]our simulated increase in arrival times from Cen A due 
to the EGMF, which averages to ~ 0.1 Myr. This is important 
for transient sources, such as gamma-ray bursts, which would 
experience at least this amount. 

We also show, in Fig. ^ the band corresponding to where 
particles with energies ^ 60 EeV produce multiple images 
of the source. We see that this runs orthogonal to our pre- 
ferred region, so that, if the presence of multiple images can 
be inferred for a given range of energies using the improved 
statistics of future data sets, a narrow range of field properties 
would be established. 

Since the appearance of multiple images can enable us to 
directly probe the properties of the local EGMF, it is tempt- 
ing to attempt to reconcile this with the specific clustering 



features seen in the Auger data; however, caution is in order 
given the present statistics. There is a total of 11 pairs of 
events in the Auger data within 5° of each other, 6 of these 
being within 18° of Cen A. Fig.[8]shows the results from sev- 
eral instances of randomly placing 13 events within a circular 
region with 18° extent in the sky. If the signal is distributed 
roughly uniformly within this region, interesting features can 
emerge from chance alone. Considering this, it is difficult to 
discern between the formation of multiple source images and 
a scattered signal at this point, although the magnetic field 
parameters yielding such distributions are not too dissimilar. 
An important question is the effect of the Milky Way's mag- 
netic field (GMF) on our results. It is thought that the GMF 
consists of a regular component with reversals in the field di- 
rection between neighboring arms of the galaxy plus a turbu- 
lent component with coherence length of ~ 0.1 kpc ( [Stanev 
1997t|Pshirkov et al. 2011). While both have field strengths 



of a few /iG, the deflection due to the turbulent component 
is considerably smaller since the regular component is coher- 
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Fig. 7. — Inferred range of extragalactic magnetic field parameters Brms 
vs. Ac (dark shaded region) that is compatible with both: (i) the average 
angular distribution of events being 8° -18° from Cen A (between the solid 
lines); (ii) the spread of events among themselves having an extent larger 
than 4° (above the dashed line). The latter condition disfavors scenarios in 
which events are shifted from the source position yet remain tightly clustered. 
These conservative requirements indicate Brms should exceed 20 nG in the 
local extragalactic environment if Cen A is the source of the Auger excess. 
As in Fig.pl the gray band corresponds to the appearance of lensing effects, 
which would allow for tighter constraints. 

ent on much larger scales. Protons with energies of 60 EeV 
are expected to be deflected by only about a degree (smaller 
than the uncertainty of UHECR detectors) and less at higher 
energies (Tinyakov & Tkac hev|2005| ). 

While these deflections could be more important for the 
regular component of the GMF, two comments are in order: 
(i) The regular component tends to produce only a coherent 
shift in t he position o f th e source (as can be se en in the sky 
maps in Stanev 1997 and'Vorob iov et al.|2009| ) and it cannot 
account for the fact that the excess events observed are scat- 
tered around Cen A; (ii) In at least some of the GMF models, 
the deflections observed from the direction of Cen A are less 
than a degree (in contrast to Galac tic Center or disk in w hich 
deflections could be much larger) ( [Takami & Sa to'20Q8 ). 

Since both types of GMF scatte ring would b e significantly 
larger for heavier nuclei (see, e.g., Giacinti et a l. 2010, 201 1 ), 
it is difficult to reconcile the observed event distribution with 
an assumption that Cen A is a source of nuclear primaries. 
This, in combination with the fact that Auger does not see 
an excess of lower-energy particles from around Cen A (as 
discussed above), suggests that if Cen A is an UHECR source, 
it is likely producing protons. We caution here that the larger 
effects of the GMF on heavy nuclei makes it non-trivial to 
directly rescale our concluded values for the nearby EGMF. 

We thus expect the GMF to have only a small effect on 
protons at the energies examined here and not to impact our 
conclusions. However, it is also possible that some as yet un- 
known magnetic field component is contributing, which could 
be even more interesting. The possible effects of the radio 
lobes or a field filling an extended halo surrounding Cen A 



will be examined elsewhere. 

7. DISCUSSIONS AND CONCLUSIONS 

Direct observations have yet to determine the strength of 
the extragalactic magnetic field near the Milky Way. While it 
is difficult to provide details of the field configuration using 
existing data, we have attempted to determine the most likely 
values under the assumption that Cen A is the cause of the 
UHECR excess coincident with its position in the sky seen by 
Auger. 

While we have only considered a Kolmogorov spectrum 
here, the overall behavior is essentially determined by the co- 
herence length of the field, so that the effects of using a dif- 
ferent power spectrum are small and do not change our basic 
conclusion that the implied local EGMF is > 20 nG. As seen 
in Fig. [T] the coherence length is degenerate with the mag- 
nitude of the field, thus for smaller coherence lengths even 
larger fields could be accommodated by the present data. 

The presence of a > 10 nG field leads to several conse- 
quences for all UHECR observations. In addition to the time 
delay discussed above, we also see that changes in the particle 
energy can drastically alter the resulting observed angular dis- 
tribution. A general feature that we find is a rapid isotropiza- 
tion with decreasing cosmic ray energy, although specific de- 
tails vary greatly between both field configurations and the 
relative positioning of Earth and Cen A. Thus inferring the 
particle spectrum by use of a fixed angular bin can be a non- 
trivial undertaking. Unfortunately, we have only a single point 
of observation, so it is quite possible that our location rela- 
tive to Cen A is such that we can observe high-energy events, 
while not receiving lower-energy cosmic rays, so that a lack 
of signal is not necessarily unexpected. This would not be 
the case if the Galactic fields were solely responsible for the 
observed dispersion, and it agrees with the lack of UHECR 
tracks aligned by energy reported by the [Pierre Auger CoUab 
oration (2012). 

Our results confirm that a simple angular projection is not 
sufficient to arrive at the nature of the local EGMF. Greater 
statistics will allow for new diagnostics, such as indepen- 
dently analyzing events from the direction of Cen A to de- 
termine their composition. If signals of magnetic lensing be- 
come evident, this would break the remaining degeneracy in 
the plane of allowed B — Ac values in Fig. [7] One expectation 
is that events in UHECR clusters caused by lensing should 
have similar energies, since otherwise a spread would result. 
Further, voids of very-low density are visible in Fig. [T] If 
these "cold spots" persist, they may yield information about 
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Fig. 8 . — Auger data from the region around Cen A are shown in the top left 
panel. This can be compared with simulated sky maps obtained by randomly 
placing 13 events within a circle of radius 18°, as shown in the remaining 
panels. This results in an average of ~ 5 — 6 pairs within 5°. 



YUKSEL, STANEV, KISTLER, & KRONBERG 



what is not producing UHECR. 

Finally, if Cen A is not a UHECR source, this would be 
important to know. It may be that the excess seen at present, 
which implies the rather large value of > 20 nG, dissipates as 
more data are obtained. The direction that future data take 
will determine what additional steps, such as a greater exam- 
ination of the effects of magnetic lensing, are warranted. 
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